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In  this  paper  we  study  the  behavior  of  a  queueing 
system  which  arises  in  the  study  of  certain  communication 
networks.  Specifically  we  study  a  queueing  phenomenon  which 
arises  with  the  SENET  network,  as  described  by  Coviello 
and  Vena  (1975)  or  Barbacci  and  Oakley  (1976) .  This  network 
allows  for  both  voice  and  data  messages  to  be  transmitted 
over  the  same  channels  by  using  a  special  type  of  integrated 
circuit  and  packet-switched  multiplexor  structure.  The  two 
classes  of  traffic  have  substantially  different  performance 
requirements.  Voice  messages  tend  to  possess  great  redundancy, 
and  hence  not  to  be  sensitive  to  channel  error  rates,  while 
data  is  very  sensitive  to  channel  error,  having  essentially 
no  redundancy.  Voice  messages  on  the  other  hand  have  critical 
timing  requirements  and  cannot  be  queued,  while  data  is 


1 


relatively  insensitive  to  timing  and  can  be  queued.  Addition¬ 
ally,  voice  messages  tend  to  be  very  long  relative  to  data 
messages  which  can  be  broken  up  into  small  packets.  These 
special  requirements  have  led  to  the  following  queueing  net¬ 
work.  A  node  of  the  network  consists  of  c  +  v  channels  or 
servers.  The  voice  messages  are  assigned  to  v  channels  and 
do  not  queue .  Thus  the  voice  messages  operate  as  a  loss 
system.  Data  messages  may  use  c  channels  exclusively  and 
any  unused  voice  channels;  however,  voice  preempts  data  using 
voice  channels.  Data  messages  are  queued  if  necessary. 

Typical  performance  measures  that  one  may  wish  to  calculate 
include  the  loss  rate  of  voice  traffic  and  the  mean  data 
queue  length. 

We  make  standard  probabilistic  assumptions.  Specifically , 
we  assume  voice  traffic  arrives  according  to  a  Poisson (A) 
process  and  each  voice  message  has  an  independent  exponential ( u) 
service  time.  Data  messages  are  assumed  to  have  independent 
exponential (n )  service  times  and  arrive  according  to  a  Poisson (6) 
process.  With  these  assumptions  voice  is  an  M/M/v/v  loss 
system,  and  data  is  an  M/M/S  system  where  S=c+v-V(t) 
with  V(t)  =  number  of  voice  messages  in  service.  The 
stochastic  process  { (X(t) ,V(t) ) ,  t  >  0}  is  Markov  with 
state  space  Z+  x  {0,1,..., v}  where  X(t)  =  data  system  sise 
at  time  t.  One  can  easily  write  the  Kolmogorov  forward 
equations  appropriate  for  this  system;  however,  these  equations 
do  not  yield  a  closed  form  solution.  To  describe  this  system 
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one  must  either  numerically  solve  the  forward  equations  or 
introduce  approximations . 

This  system  has  been  studied  previously  by  a  number 
of  researchers  including  Halfin  and  Segal  (1972) ,  Halfin  (1972) , 
Fischer  and  Harris  (1976),  Bhat  and  Fischer  (1976),  Fischer 
(1977) ,  Chang  (1977) ,  and  Gaver  and  Lehoczky  (1979a, b) .  The 
last  two  papers  introduce  a  "fluid  flow"  and  a  diffusion 
approximation  and  derive  explicit  formulas  for  data  queue 
behavior.  These  papers  focus  on  the  important  case  in  which 
Pd  =  6 /n  >  c.  In  such  a  situation  the  data  messages  must 
have  access  to  voice  channels  for  the  system  to  be  stable. 

4 

Furthermore,  it  was  assumed  that  n/p  was  large,  say  10  . 

Under  these  circumstances  the  data  flow  could  be  treated 
deterministically.  Suppose  we  define  pv  =  X/p  and 
q  =  (Py/vl)/[j_0  P^/j! ,  the  Erlang  B  blocking  probability. 

The  total  traffic  intensity  on  the  c  +  v  channels  is  given 
by  pd  +  pv(l-q),  or  we  could  define  p  =  (od  +  pv(l-q) ) / (c+v) . 
A  heavy  traffic  approximation  can  be  derived  for  this  case 
p  f  1.  Such  an  approximation  was  derived  in  Gaver  and 
Lehoczky  (1979b)  assuming  n/p  was  large;  a  Wiener  process 
with  reflecting  boundary  was  found  appropriate.  In  this 
paper  we  derive  a  heavy  traffic  approximation  for  the  system 
without  the  fluid  flow  assumption  that  n/p  is  large. 

The  methodology  is  drawn  heavily  from  the  approach  of  Burman 
(1979) .  In  this  approach  one  characterizes  a  Markov  process 


by  its  infinitesimal  generator.  One  next  suitably  normalises 
the  process  so  that  the  generator  converges  to  a  limiting 
infinitesimal  generator  (in  this  case  to  that  of  a  reflected 
Brownian  motion) .  This  convergence  allows  the  conclusion  that 
the  finite  dimensional  distributions  of  the  normalized  Markov 
process  converge.  The  diffusion  approximation  consists  of 
treating  the  actual  process  through  its  limiting  behavior. 

The  details  are  somewhat  complicated  by  the  presence  of  a 
boundary . 


2. 


'y* 


Let  { (X(t) ,V(t) ) ,  t  _>  0}  be  a  bivariate  Markov 
process  with  state  space  S  »  Z+  x  (0,1, . . . # vJ  •  Here 
{V(t) ,  t  >  0}  is  marginally  an  M/M/v/v  loss  system  with 
arrival  rate  A  and  service  rate  y.  Conditional  on  V(t) , 
{X(t),  t  >  0}  is  an  M/M/ (c  +  v  -  V(t))  queueing  system  with 
arrival  rate  6  and  service  rate  n .  We  say  that  the  V 
process  subordinates  the  X  process.  We  let 


the  infinitesimal  generator  of  the  V  process. 

The  generator  of  the  (X,V)  process  is  given  by 


Qf(x,k)  +  5(f(x+l,k)  -  f(x,k)) 

+  n(c+v-k) (f (x-l,k)-f (x,k) ) 


Af{x,k)  = 


if  x  >  c+v-k 


Qf(x,k)  +  6(f(x+l,k)  -  f  (x,k) ) 

+  nx(f(x-l,k)  -  f(x,k)) 

if  x  <  c+v-k 


(2.2) 
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for  f:S  ■*  R  continuous  where 


Qf(x,k)  -  pvf(x,k+l)  -  (k+pv)  f(x,k)  +  kf(x,k-l) 

v  >  k  >  C 


(2.3) 


and  f (x,-l)  ■  f(x,v+l)  =  0.  Clearly  Qf(x)  »  0,  that  is  Q 
annihilates  functions  of  x  alone.  We  next  normalize  the 
(X, V)  process  by  defining  Xft(t)  *  X(nt)//n  and  Vn(t)  =  V(nt)  . 
One  can  calculate  the  generator  of  the  Markov  process 
{  ( Xn ( t ) ,  Vn(t)),  t  2  0}  having  state  space 
Sn  =  (0,  2//n , .  .  .  }  x  {  0 ,1 ,  . . . ,  v)  to  be 


Aflf (.x,k) 


nQf (x,k) 


nQf (x,k) 


+  <5  n(f(x  +  l//rf,k)  -  f(x,k)) 

+  nn(f(x  -  1/Vn,k)  -  f(x,k)) 

if  x  >  C  *  V  k  (2 

/n 

+  5n(f(x  +  l//”n,k)  -  f(x,k)) 

+  nn  /nx  ( f ( x  +  l//n,k)  -  f(x,k)) 

if  x  =  0,  l//n,  ...,  (ctv-  k) //H  . 


We  assume  f(x,k)  has  three  bounded  derivatives  in  x 
for  each  fixed  k.  With  this  assumption  one  can  expand  terms 
in  (2.4)  in  a  Taylor  series  and  rewrite  as 
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r 


ARf (x, k) 


- 


nQf(x,k)  +  n1‘//2fx(x,k)  (5 -n  (c+v-k) ) 

+  ^  f^fXfk)  (6 +n  (c+v-k) ) 

♦  (Xn-1'2) 
if  x  (c+v- k)//n 

nQf(x,k)  +  n^2fx  (x,k)  (6 -r)/nx) 

+  ^  f„„(x,k)  (6+n  /nx)  +  o(n“1/,2> 


xx 


(2.5) 


if  x  =  0,  l//rT,  ...  ,  (c+v-k)//n 


with  fx(x,k)  =  f(x,k)  and  fxx(x,k)  =  — 2-jf(x,k). 

3x 

We  ultimately  wish  to  prove  that  the  finite  dimensional 

distributions  of  (Xn(t),  t  )  0)  converge  to  those  of  a 

Wiener  process  with  reflecting  barrier  at  the  origin.  This 

can  be  restated  in  terms  of  semi -groups .  We  let  (t",  t  ^  0} 

be  the  semi-group  of  operators  associated  with 

{ (X  (t) ,  V  (t)),  t  >  0}  and  {T~,  t  >  0}  be  that  associated 
n  n  —  t  — 

with  a  Wiener  process  having  reflecting  barrier  at  0.  Let 
g  be  a  continuous  function  g:R'  -*■  R'  .  Knowledge  of  the 
semi-group  is  equivalent  to  knowledge  of  the  transition 
functions  by  taking  a  sequence  of  g's  which  approximate 
indicator  functions.  We  wish  to  prove  |t£  g(x,k)  -  T~  g(x) I  ♦ 
as  n  -*■  «  for  all  (x,k) .  Here  g(x,k)  =  E(g(Xt) 1x^(0)  = 

Vn(0)  =  k) .  The  presence  of  the  variable  k  prevents  this 


x 


from  being  done  directly.  The  method  we  use  is  to  construct 
a  convenient  sequence  of  functions  <gn>n_^  which  converge 
in  some  sense  to  g.  We  write 


I1?9  -  T;g»n  <  !)<T~g)n  -  T~g|ln  +  lITjg  -  T^gn!ln 

+  i|T?9n  -  <V>nl!n 


(2.6) 


where 


I  refers  to  the  sup  norm  over  S  .  Both  t"  and 
n  *  n  u 


Tfc  are  contraction  semigroups . 

<(T*g)  >°°  .  is  the  sequence  of  functions  constructed 
t  n  n=i 

from  T^g.  Our  goal  is  to  show  that  each  of  the  three  terms 

on  the  right  side  of  (2.6)  converges  to  0.  The  first  and 

second  terms  can  be  handled  similarly.  For  any  function  g, 

we  must  guarantee  that  the  constructed  <9n>n=i  sequence 

satisfies  ||g„  -  gIL  -*■  0.  It  will  follow  that  l!(T~g)  -  T~g|l  H 
n  n  t  nt  n 

Moreover,  since  (T^,  t  >_  0}  is  a  contraction  semi-group 

liT^g  -  TVg  l|  <  ||g  -  g  II  which  also  converges  to  0.  The 
t  t  n  n  —  n  n 

sequence  <gn>~_-^  he  chosen  in  such  a  way  that  the  third 

term  converges  to  0 . 

We  focus  on  a  convergence  determining  class  of  functions 
g,  those  which  are  bounded  and  have  three  bounded  derivatives. 
For  such  a  function  g(x)  we  define 


gn(x,k)  =  g(x)  +  — u(x,k)  +  i  v(x,k) 
1  /"n 


(2.7) 
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where  u  and  v  have  two  bounded  derivatives  in  x  for 

each  fixed  k.  The  functions  u  and  v  will  be  determined 

explicitly  later  and  are  chosen  to  control  the  third  term  in 

(2.6).  Clearly  when  g  is  defined  by  (2.7),  ||g  -g||  =  o(n-1^2) 

n  n  n 

and  therefore  converges  to  0  as  required. 

One  can  apply  the  generator  An  to  gn  to  derive 


Angn(x,k) 


jnQq(x)  +  n1//2  [Qu(x,k)  +  g'  (x)  (6  -n  (c+v-k) )  ] 
+  [Qv (x,k)  +  ux(x,k)  (5 -n  (c+v-k) ) 

+  \  g"  (x)  (6  -n  (c+v-k) )  ]  +  0(n"1/2' 
if  x  >  (c+v-k)  //n 


~  \ 


(nQg  *  *  +  n1//2  [Qu(x,k)  +  g’  (x)  (6-  nn1//2x'  ] 

+  [Qv(x,k)  +  u^  (x,k)  (5  -nn1^2x)  +  ^  g"  (x)  (S  +nn1/^2x)  ] 


+  0(n-1/2) 


if  x  =  0,  l//n,...,  (c+v-k) //n 


g 

where  ux(x,k)  =  u(x,k) .  Recall  that  Q  annihilates  functions 

of  x  alone,  thus  nQg(x)  =0.  We  want  to  have 

converge  to  a  finite  limit  and  to  have  that  limit  be  inde- 

1/2 

pendent  of  k.  For  this  to  occur,  the  n  '  term  must  be 
controlled  and  the  functions  u  and  v  must  be  chosen  in 
such  a  way  as  to  eliminate  the  variable  k. 
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1/2 

The  n  '  coefficient  in  (2.8)  can  be  rewritten  by 
adding  and  subtracting 
v 

l  irvg*  (x)  (6  -X]  (c+v-k) )  =  -n  (c+v)  (l-p)g' (x)  . 
k-0  * 

We  next  pick  u(x,k)  to  be  a  solution  of 

Ou(x,k)  =  - (g‘ (x)  n (Pd~ (c+v-k) )  +  g' (x)  n(c+v)(l-p)) 

=  -g' (x)  n(k-pv(l-q))  (2.9) 

When  u(x,k)  is  any  solution  of  (2.9),  the  coefficient  of 
1/2 

the  n  '  term  in  (2.8)  becomes 

-g'(x)  n  (c+v)  (1-p)  if  x  >  -+-V-^ 

/n 

g'(x)  n ( (C+v) p  -  n1/2x  -  k)  if  0  <  x  <  cf?~l 

/n 

Equation  (2.9)  can  be  solved  explicitly.  Define 
ak  =  -g' (x)  n(k-Pv(l-q) ) /y,  so  that  (2.9)  can  be  written  as 

-Pv(u(x,k) -u{x,k-l) )  —  (k-1 )  (u(x,k-l) -u(x,k-2) )  =  k  *  l»...»v 

-v(u(x,v)-u(x, v-1) )  *  av  (2.10) 
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Equation  (2.10)  has  a  solution  since  JkmQ 
where  <n.  >Y  n  is  the  stationary  distribution  associated 

K  K=U 

with  Q,  or  7rk  *  (p^/ki )/^Y[-q  P^/i'.).  The  solution  is 
given  by 


u(x,k)  -  u(x,k-l)  = 


k-1 

I  u  .a. 
i=0  1  1 

0  7T,  i 

Fv  k-1 


-g‘ (x)nTk_1 


VfVk-l 


where 


Tv  *  l  n.{l-p  (l-q)>  and  T  =  0. 
*  i=0  1  v  v 


Clearly 


u(x, k)  »  u(x,0)  - 


j  Ii=I  ,  i  <  k  < 

UP_.  L  71  ,  —  - 


v  (2.11) 


v  t»i  i'1 


where  u(x,0)  is  arbitrary.  We  let  u(x,0)  *  ^  g* (x)  so 


K  T  • 

u(x,k)  *  g’  (x)  \  |  -  -g-  l  ^  |  .  0  <  k  <  v  (2.12) 

MMv  i*l  i-] 


For  the  choice  of  u  specified  by  (2.12)  we  next  wish 


to  insure  that  the  limiting  generator  is  independent  of  the 
variable  k.  The  function  v  is  chosen  to  eliminate  the 
dependence  on  k.  The  0(1)  term  of  (2.8)  is  given,  for 


Qv(x,k)  +  g"(x)  i  -  l  -p-i  (6  -n  (c+v-k) ) 

I  1  wpv  i«l  *i-l  J 

+  j  gH  (x)  (S  +n  (c+v-k) ) 

»  Qv(x,k)  +  H (x,k) . 

Let  H(x)  =*  IfcaQ  7TkH(x,k)  and  consider  Qv(x,k)  +  (H  (x,k)  -H (x) )  +H (x)  . 
We  now  let  v(x,k)  be  any  solution  of 

Qv(x,k)  =  - (H (x,k)  -  H ( x) )  .  (2.13) 


Equation  (2.13)  has  a  one-parameter  family  of  solutions, 
since  £k=Q  ^k(H(x,k)  -  H(x))  *  0.  When  v(x,k)  is  chosen 
to  be  any  solution  of  (2.13),  the  0(1)  term  of  (2.3), 

for  x  (c+v-k)  //n,  will  become  H(x)  and  will  therefore 
be  independent  of  k.  It  remains  to  calculate  H(x) . 


H(x)  =  g" (x) 


[!•■  |(s  -  *  j,  v;) 

')] 


+  2  (5  +n  (c+v-k) ) 


g"  (x) 


s'  ^  Jo  j,  ft'1 


i=l  i-1 


9"(x)  [“4  Jo  "’‘<k-pv(l-q”  Jx 


k  t. 


+  n (c+v) (1-p)  y  y  *i-l 

UP„  k“0  k  ir4_. 


i*l  "i-l 


(2.14) 
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The  second  term  can  be  rewritten  by  interchanging  the  order 
of  svanmation.  The  third  term  is  O(l-p).  We  find 


H(x)  =  g" (x) 


g"  (x) 


2  v— 1  T.  v  I 

6--°-  l  —  l  *k(k-p  (1-q))  +  O(l-p) 

L  upv  i=0  »i  k=i+l  k  v  J 

2  v-1  T .  I 

*7  (Tv'V  +0,1-°,J 


with  Tv  =  0  or 


H(x)  =  g"  (x)  t\ 


v-1  T? 

pd  +  £§“  .L  FT  +  0(1“p) 

v  1=0  i 


(2. IS) 


For  the  functions  u  and  v  specified  by  (2.12) 
and  (2.13),  equation  (2.8)  can  be  rewritten  as 


Angn(x,k) 


-n1^2  (i-p )  (o+v)  ng’  (x) 
T  l  V_1  Ti 

+  Td  +  v—  Jo  7 


+  o(l-p) 


g" (x)  +  o(n~1/2) 


for  x  >  (c+v-k)//n 


(2.16) 


n1//2  t\  ( (c+v)  p  -n1//2x-k]  g' (x) 

v-1  T? 


+  n9"(x) 


pd  +  pj"  l  TT  +  0(1‘p) 

v  1=0  i 


-  (c+v 


-k-n1'2*)  -H-  I  1*1 

ppv  iio  "i  J 


+  Otn-1'2) 


fOr  x  <  (c+v-k)//n 
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4 

i 

x 


We  now  introduce  the  "heavy  traffic  approximation." 

In  order  for  the  generator  to  converge  to  a  limiting 

-1/2 

generator  we  must  have  1-p  *  0(n  '  ) .  Specifically,  we 
assume  p  ■  p  ■  1  -  (0//n)  for  some  6  _>  0 .  In  this  case, 
n1/^2(l-p)  *  0,  and  (2.16)  becomes 


-0n (c+v) g' (x)  +  n 


r  „  v-i  t2i 

h +  %  Jo 


g" (x)  +  0(n"1/2) 


for  x  >  (c+v-k)//n 


A  g  (x,k) 
n  n 


(2.17) 


n  I  (c+v)  p-n1^2x-k]n1/^2g'  (x) 


+  ng"  (x)  N  p  .  + 


-  (c+v-k-n1/^2x) 


k  T.  1 

y  _± 

PV  do  \  J 


for  x  £  (c+v-k)//n 


We  now  define  a  limiting  generator  Am  with  domain 
consisting  of  all  functions  g  having  three  bounded  derivatives 
and  g' (0)  *  0.  Let 


A  g(x)  =  -0rj  (c+v)  g'  (x)  +  n 


r  v-i  t2  i 

Pa  +  —  l  —  g"  (x)  , 
L  d  Upv  i*0  \] 


x  >  0 


(2.18 
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■  V'.  ev’  i Ti W  v ■» 


Aw  is  the  generator  of  a  Markov  process  which 
corresponds  to  a  Wiener  process  with  drift  -8n(c+v),  scale 


2ri  pd  +  *r 


v-l  T2  1 

l  r-  I- 

i»0  *i  J 


-1/2 

and  a  reflecting  barrier  at  0.  The  0(n  )  terns  involve 

the  first  three  derivatives  of  g  which  are  bounded.  It  is 
clear  from  a  direct  comparison  of  (2.16)  and  (2.18)  that 
| A  g  (x,k)  -  A  g|  +  0  as  n  -*■  «  for  all  x  >  0  and  k 
arbitrary.  In  addition,  g1 (0)  =0  is  necessary  for  the 
generator  to  converge  at  x  -  0.  unfortunately  even  assuming 
g'  (0)  =  0, 

2  k  T2 

lAng(0,k)  -  A^gtO)  |  -  (c+v-k)  ~  g"(0)  as  n  -  - 


rather  than  to  0.  One  needs  a  special  argument  to  handle 
this  lack  of  convergence  at  the  boundary. 

We  set  out  to  prove  the  third  term  in  (2.6)  converges 
to  0.  Standard  semi-group  results  (see  Burman,  1979,  p.  33) 


give 


(T~g)  -  T"g  =  /  Tj  _  (  (A  w)  -  A  w  )  dS 

t  n  t  n  q  t-s  00  n  n  n 


(2.19) 


where  w  =  w(t,x)  =  Tfcg(x) .  Recall  that  wR  *  w  + (l//n) u+(l/n) v 
with  u  and  v  defined  by  (2.12)  and  (2.13)  with  g  replaced 
by  w.  It  follows  that 


15 


»*!«»-<*>>  A 

'  »/  <.s  <<*.“>  n-d.W+A.W-An»n)dSlin 

<  /'  »T?.sIIA.«)n-A.w)»n<3S  +  II ^.sU.W-Anwn)<JS|ln 

t  t 

<  /  II  (Aw)  -Aw  |t  dS  +  ||/  TV  (Aw-Aw)  dS||  • 

—  q  «>  n  “>  n  q  t— a  °°  n  n  n 

-1/2 

The  first  term  is  clearly  0(n  ) .  It  remains  to 

-1/2 

show  that  the  second  is  0(n  )  as  well.  We  have  shown 

I A  w  -  A  w  |  »  0(n-1/2)  except  at  the  boundary  where  it 

is  0(1) .  We  split  the  integral  into  two  parts,  for  one 

of  which  the  process  is  away  from  the  boundary,  and  for  the 

other,  near  the  boundary.  The  integral  away  from  the  boundary 

-1/2 

has  an  integrand  which  is  0(n  '  ) .  The  integral  near  the 

-1/2 

boundary  is  also  0(n  '  )  since  under  a  heavy  traffic 

assumption  the  process  is  rarely  near  the  boundary.  The 

details  are  merely  summarized  here;  they  are  based  on  the 

ideas  of  Burman  (1979). 

Let  I  be  the  indicator  function  of 
on 

To.  2^1 

L  /n 

and  be  the  indicator  of 
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We  have 


tc+y-k 
✓n 


-) 


1/  T?.s(A.w  -  Vn)ds"n 


<  II  f  T^.s(A.w  -  V„)Iln<iSlln  *  ll/  ^.S(».W  -  AnwnII0ndS!ln 

<  ll/1  <»,«  -  Vn11!,®*,  +  «»-»  -  H  T?-sIOndSllr,  • 


The  first  term  is  0(n  ,  since  |A  w  -  Awl  *  0(n  J'/*)  off 

00  n  n 

the  boundary.  The  factor  II A  w  -  A  w  ||  *  0(1),  thus  it  remains 

00  n  n 

to  show  that 


.“1/2, 


11  /  Tt-sI0nds^n  -  0("'1/2» 


This  gives  the  total  time  in  (0,t)  spent  near  the  boundary. 
We  bound 


"/  T?-SI0ndSlln 


by  first  introducing  a  function  h(x)  not  in  the  domain  of 

Aoo.  We  let  h(x)  have  bounded  support,  be  infinitely 

differentiable  and  be  given  by  h(x)  =  x  for  x  near  0. 

One  can  construct  hn(x)  using  (2.7)  and  apply  An  to 

h  to  find 
n 


17 


0(1) 


if  x  > 


c+v-k 

/n 


A_h 
n  n 


(2.20) 


n1/2n((c+v)p-n1/2x-k)  +0(1)  if  x  < 

/n 


One  has 


T?h  -  h  *  /  tJa  h  dS 
t  n  n  S  n  n 


*  l  TSAnhnIlndS  +  [  feV1' 


It  follows  that 

“Z1  ❖nV^'n  i  KV  Vn  +  ’sAnhnIlndS»n 


<  2t|hnlln  +  0(1)  . 


We  have  shown  I!/  TSAnhnIondA|ln  to  be  bounded  in  n*  An 
application  of  (2.2)  shows 

H/1  TsVn'on^V  *  "l/2"  I  «»v>  p  -  "l/2x  -  k  +  0111  1 11/  Vo.>d,,lln 


is  bounded  in  n.  It  follows  that  \\f  Tg1ondS^n*'  °^n  • 
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This  finally  concludes  the  argument  which  shows 

llTt9n  "  (Ttg> n^n  *  °<n‘1/2^  hence  by  (2*6>  HT?9  -  V||n  «  0{n~1/2) 
We  have  thus  shown  that  the  finite-dimensional  distri¬ 
butions  of  the  (Xn(t) ,  Vn(t))  process  converge  to  those 
of  a  Wiener  process  with  drift  -6n(c+v)  scale 


n 


v-l 

l 


i=0 


9 


and  reflection  at  0.  The  diffusion  approximation  treats 
Xn(t)  as  though  it  were  such  a  Wiener  process.  For  instance, 
the  limiting  wiener  process  has  a  stationary  exponential  dis¬ 
tribution  with  parameter 


0  (c+v) 
v=T"“ 


pd  + 


-2—  I  (TV*-) 

wpv  i«0  1  1 


This  is  a  distribution  for  X(nt)//n  and  suggests  X(t) 
will  have  a  steady  state  distribution  given  approximately 
by  an  exponential  with  parameter 


(c+v) (1 


-/(•■•  ^2  ?)• 


The  steady  state  mean  data  queue  length  would  then  be 
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*; 


the  results  of  Gaver  and  Lehoczky  (1979b)  since  the  variability 
in  the  data  queue  is  now  included.  When  n/u  is  large,  the 
second  term  dominates,  and  the  fluid  flow  approximation  is 
satisfactory. 

The  Wiener  process  approximation  for  the  X(t)  process 
provides  a  method  for  studying  the  dynamics  of  that  process. 

For  instance,  suppose  the  data  queue  were  at  level  x  at 
time  t  where  x  is  large.  One  might  wish  to  study  the  time 
that  elapses  until  the  queue  becomes  empty.  This  is  essentially 
the  duration  of  the  busy  period  under  heavy  traffic  and 
corresponds  to  a  first-passage  time  for  a  Wiener  process.  Let 
us  denote  it  by  T  .  Straightforward  martingale  arguments 
provide  for  its  transform 


-sT 

E  (e  x)  =  exp 


(£)  -  (£> 

o  0 


(~)  +  2s 


(2.22) 


where 


m  =  9(c+v)n  ~  n^2  (l-p)  (c+v)n 


a2 

T 


,  „  n-1  t. 

n  p  +  _H_  y 

1  d  ypv  i-0  *i 


It  is  also  easy  to  find  the  mean  first-passage  time 


E(TX)  =  x/m 


(2.23) 
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One  might  also  be  interested  in  the  area  beneath  the 
sample  path  until  emptiness  occurs,  since  this  area  represents 
the  total  time  waited  by  all  data  customers  involved  in  the 
busy  period.  If  Ax  represents  this  area,  simple  backward 
equation  arguments  give 


E(Ax> 


where 


m  and 


2 

o 


are  given  in  (2.22). 


(2.24) 
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